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Abstract 

We calculate the expected flux of 7-ray and radio emission from the LMC due to 
neutralino annihilation. Using rotation curve data to probe the density profile and 
assuming a minimum disk, we describe the dark matter halo of the LMC using 
models predicted by N-body simulations. We consider a range of density profiles in- 
cluding the NFW profile, a modified NFW profile proposed by Hayashi et al. (2003) 
to account for the effects of tidal stripping, and an isothermal sphere with a core. We 
find that the 7-ray flux expected from these models may be detectable by GLAST 
for a significant part of the neutralino parameter space. The prospects for exist- 
ing and upcoming Atmospheric Cherenkov Telescopes (ACTs) are less optimistic, 
as unrealistically long exposures are required for detection. However, the effects of 
adiabatic compression due to the baryonic component may improve the chances for 
detection by ACTs. The maximum flux we predict is well below EGRET's mea- 
surements and thus EGRET does not constrain the parameter space. The expected 
synchrotron emission generally lies below the observed radio emission from the LMC 
in the frequency range of 19.7 to 8550 MHz. As long as (o~v) < 2 x 10 _26 cm 3 s _1 
for a neutralino mass of 50 GeV, the observed radio emission is not primarily due 
to neutralinos and is consistent with the assumption that the main source is cos- 
mic rays. We find that the predicted fluxes, obtained by integrating over the entire 
LMC, are not very strongly dependent on the inner slope of the halo profile, varying 
by less than an order of magnitude for the range of profiles we considered. 
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1 Introduction 



Decades of observational evidence indicate that galaxies are surrounded by 
massive dark matter halos. It is also well-established that the baryon fraction 
of the universe is too small to account for the measured dark matter density. 
Most recent studies suggest that ~ 23% of the energy density of the universe 
is in the form of non-baryonic cold dark matter (CDM) (see, e.g., [1]). 

Of the many candidate dark matter particles that have been proposed, the 
most popular is the Lightest Supersymmetric Particle (LSP). Strongly moti- 
vated by theoretical considerations, the LSP is stable in R-parity conserving 
models and has the appropriate relic density for masses m LS p ~ 100 GeV 
(see, e.g., [2]). In Minimal Supersymmetric extensions of the Standard Model 
(MSSM), the LSP is most likely the lightest neutralino, x [3]- Under certain 
assumptions, WMAP results favor m x < 500 GeV [4], while accelerator ex- 
periments constrain m x ^ 50 GeV [5]. 

Neutralinos may be detected directly in low background experiments or indi- 
rectly by observations of their annihilation products. Products of XX annihila- 
tion rapidly decay into neutrinos and 7-rays, as well as electrons and positrons 
which emit synchrotron radiation in local magnetic fields. The rate of XX an_ 
nihilation depends on the square of the neutralino number density. As such, 
the largest fluxes are expected to come from the densest regions. The galac- 
tic center is thus an obvious candidate for its density and proximity [6,7,8], 
but the numerous associated background signals might make it difficult to 
disentangle the neutralino annihilation flux. 

High resolution N-body simulations of CDM halos find an abundance of sub- 
structure in galaxy-size halos. These studies predict that satellites of our 
galaxy, such as dwarf spheroidal galaxies, are embedded in dark matter halos of 
their own. In addition, these simulations find numerous dark "clumps". The 
detection of neutralino annihilation products from these substructures may 
provide information about the clump density profiles as well as neutralino pa- 
rameters. Previous studies have considered the annihilation signal from dwarf 
spheroidal galaxies [9,10,11] and dark matter clumps [12,13,14,15,16,17,18]. 
The detectability due to neutralino annihilation of other systems such as the 
LMC [19], M31 [20], and M87 [9], the globular cluster Palomar 13 [21], and 
galaxy clusters [22] have also been studied. 

At a distance of 50.1 kpc [23], the LMC is an obvious choice for neutralino de- 
tection. There are several LMC mass estimates in the literature (e.g., [24,25,26,27,28,23]), 
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and most agree within error bars. According to the most recent study, the mass 
of the LMC enclosed in a radius of 8.9 kpc is M LMC ~ (8.7 ± 4.3) x 1O 9 M . 
To study neutralino annihilation it is important that the LMC be dark matter 
dominated. Van der Marel et al. [23] found that more than half of the LMC 
mass is in a dark halo. In addition, Sofue [27] argued that the LMC must 
have a dark and compact bulge with an anomalously high mass-to-light ratio 
M/L ~ 20 - 50 M & /L & . 

The LMC has been observed in a wide range of frequencies, providing a wealth 
of data to draw upon. Past estimates of the LMC neutralino signal [19] focused 
solely on 7-rays and did not fit modern density profiles to the LMC halo. New 
rotation curve data and recent theoretical models of tidally stripped dark 
matter sub-halos motivate a new estimate of the LMC neutralino flux in 7- 
rays and in synchrotron emission. In addition, next-generation Atmospheric 
Cherenkov Telescopes (ACTs) and the upcoming GLAST satellite make this 
study of the LMC even more timely. 

The flux of both 7-ray and synchrotron emission from neutralino annihilation 
in the LMC depends on the square of the density profile of the LMC dark 
halo. Thus, we begin the next section with a discussion of the structure of the 
LMC and describe fits of different dark matter profiles to the LMC rotation 
curve data. In §111 we calculate the 7-ray flux and investigate the detectability 
prospects for ACTs and GLAST. In §IV we calculate the synchrotron flux and 
compare it to observations. We conclude in §V. 



2 Modeling the LMC dark matter halo 

Rotation curve data provides a dynamical measure of the mass distribution 
via 

t/ ( \ \ GM{r) 

V rot{r) = \j , (1) 

where V rot (r) is the rotation velocity at radius r and M(r) is the total mass 
interior to r, M(r) = J r p{r )d 3 r . Given an observed rotation curve, we can 
study how well a particular model of the density distribution fits the observa- 
tions. For the LMC rotation curve data, we use HI data from Kim et al. [26] 
that spans the range from 0.05 to 3.09 kpc with a velocity resolution of ~ 1.65 
km/s, and carbon star data from Alves and Nelson [28] that covers a radius 
range of 4 to 8.2 kpc. We assume 8% error bars for the HI data given in [26]. 

Numerous N-body simulations have focused on understanding the structure 
of dark matter halos. Navarro, Frenk and White (NFW) [29,30] performed 
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N-body simulations of CDM halos and found that for a large range of masses 
the density profiles were well-described by the simple formula 



p(r) 



Po 



(2) 



(r/r s ) (l + r/r s ) 



2 ' 



with p and r s a characteristic density and scale radius, respectively. The 
NFW profile can be alternatively characterized by r s and the concentration 
parameter c which is related to p Q and is defined as the ratio r v /r s where r v 
is the virial radius. The virial radius is defined as the radius within which a 
certain virial overdensity is reached, typically 180 times the average density of 
the universe. Behaving as r^ 1 at the center, the NFW profile is shallower than 
both the singular isothermal sphere and the Moore et al. [31] density profiles 
which behave at the center as r~ 2 and r~ L5 , respectively. Given the apparent 
discrepancy between the predictions of N-body simulations and observations 
(see, e.g., [32]), with the latter favoring shallow rather than very steep profiles 
at the center, we primarily focus on the NFW profile. This is a conservative 
assumption for CDM profiles since steeper density profiles lead to larger fluxes 
from neutralino annihilation. In particular, the effect of the baryonic compo- 
nent on the central density profile may increase the central density by more 
than an order of magnitude (see, e.g., [33]). 

Hayashi et al. [34] performed N-body simulations to investigate how the den- 
sity profile of an NFW halo changes due to tidal stripping. They found that 
the resulting profile can be described by a simple modification to the NFW 
profile: 



The factor f t is a dimensionless measure of the reduction in central density 
while the "effective" tidal radius r te describes the outer cutoff due to tides. The 
effect of tidal stripping on the initial NFW profile is to lower the characteristic 
density, to make it steeper in the outer region, and to introduce an outer cutoff. 
Since there is evidence that the LMC may be undergoing tidal disruption, we 
also model its halo using this profile. 

As observations find that more than half of the mass within ~ 9 kpc of the 
LMC is dark [23], we assume that the disk is minimal (i.e., that the gravi- 
tational potential of the LMC is entirely due to dark matter) and search for 
the best fit NFW and Hayashi et al. profiles. For the NFW profile we find 
r s = 9.16 kpc and p Q = 1.66 x 10 7 M Q /kpc 3 , corresponding to a concentration 
index c ~ 11.2, and a x 2 P er degree of freedom of 1.01. 



p(r) 



1 + (r/rte) 



3 Pnfw " 



(3) 
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For the purposes of fitting the Hayashi et al. profile, we absorb the factor f t 
into the NFW parameter p and call the resulting parameter p . To distinguish 
between different models with almost identical \ 2 P er degree of freedom values, 
we first fix the scale radius at the NFW best-fit value, r s = 9.16 kpc. In this 
case, we find r te = 6.12 kpc, p Q = 8.18 x 10 6 M Q /kpc 3 with a \ 2 P er degree of 
freedom equal to 0.91. 

Recent N-body simulations ([35], see also [34,36]) suggest that the process of 
tidal disruption of an NFW halo of mass comparable to the LMC may lead to 
a decrease in the scale radius of up to about 30%. The density profile of the 
LMC at the present is well-fitted by r s = 9.16 kpc, so we can assume that this 
value represents a decrease of not more than 30% from the original value. Thus 
we can extrapolate that the scale radius of the original NFW halo that gave 
rise to the LMC was not more than ~ 13.1 kpc (~ 9.16/0.7 kpc ~ 1.4 x 9.16 
kpc). To understand the effect that a larger scale radius has on the Hayashi 
et al. density profile and the resulting flux, we again fit the rotation curve 
data using a fixed r s = 13.1 kpc. The resulting best-fit Hayashi et al. profile 
parameters are given by r te = 4.97 kpc, p Q = 5.46 x lO 6 M /kpc 3 with a x 2 
per degree of freedom equal to 0.81. 

All three fits to the rotation curve data are shown in Fig. 1. While all the 
fits are acceptable, the Hayashi et al. profiles better fit the data at the outer 
radii. The masses predicted by these models are in reasonable agreement with 
estimates in the literature - the Hayashi et al. gives a mass of ~ 9 x 1O 9 M 
within 8.9 kpc, while the NFW fit gives ~ 3 x 1O 1O M within 8.9 kpc, slightly 
above the observationally determined (8.7 ± 4.3) x 1O 9 M within 8.9 kpc 
by [23]. 

The relatively small tidal radii of the Hayashi et al. fits deserve comment. 
Observations typically indicate the tidal radius for the LMC r t is greater than 
about 10 kpc ([23] and references therein). However, different observational 
methods for determining this quantity may give significantly different values. 
While a precise determination of the tidal radius would require knowing the 
distance at which regular isopleths of a certain class of stars end, more practical 
estimators involve the Roche limit (valid for point masses), or the Jacobi limit 
(where the centrifugal force on the satellite is taken into account). For instance, 
the tidal radius determined via the Roche limit (r t R = [Mlmc /2Mmvk) 1//3 ^) 
relies upon approximate mass estimates for both the Milky Way halo (Mmw) 
and the LMC. These methods do not directly correspond to the Hayashi et 
al. profile parameter r te , which is an effective tidal radius. Careful inspection 
of the enclosed mass as a function of radius for tidally stripped objects in 
the Hayashi et al. study reveals that the mass does not necessarily reach a 
plateau at r = r te , but continues to increase with radius (e.g., up to 2 times 
r te , which is closer to the observational r t ). Consequently the parameter r te 
does not necessarily correspond to the observationally determined tidal radius. 
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Fig. 1. Rotation curve of the LMC: square points from HI data [26] and triangle 
points from carbon star data [28]. The crossed curve is the NFW fit. The dashed 
curve is the Hayashi et al. fit obtained by fixing the scale radius to the best fit 
value found for the NFW profile. The solid line corresponds to the Hayashi et al. 
fit obtained by fixing the scale radius at r s ~ lAr s ^FW to account for the possible 
evolution of r s (see text for details). 

Lastly, even if we force r te = 15 kpc, we find that the Hayashi et al. fit becomes 
essentially the same as the NFW best fit (the Hayashi et al. best fit r s = 9.35 
kpc while the NFW best fit r s = 9.16 kpc). This is not surprising given that 
with such a large r te the difference between the Hayashi et al. and the NFW 
profiles arises at radii much larger than that for which there is rotation curve 
data available. 

For completeness, we consider three additional profiles. We examine two re- 
cently proposed profiles: the Stoehr et al. [37] profile which was found to 
accurately fit rotation velocity curves of satellite halos in some simulations, 
and the more centrally concentrated Moore et al. [31] profile. We also look at 
a shallower profile - an isothermal sphere with a core as derived in [28] by 
fitting the same rotation curve data we use here. Both the Moore et al. [31] 
and the Stoehr et al. [37] profiles seem less favored by the rotation curve data 
- they give x 2 values per degree of freedom 5.49 and 4.24 respectively. The 
isothermal sphere with a core represents a less ideal scenario. Instead of mak- 
ing the assumption of a minimum disk, the stellar disk and gas contributions 
are included in the model by Alves and Nelson [28]. In this case, the halo 
profile is given by 
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with p = lO 8 M /kpc 3 and core radius a = 1 kpc. 



3 The 7-ray Emission 



Neutralino annihilation proceeds through a number of channels, many of which 
produce 7-rays in the final state (see e.g., [8]). Here we consider only the con- 
tinuum 7-ray emission due to the decay of neutral pions produced in hadronic 
jets. 



^7T°^7 + 7. (5) 

Given the uncertainties in the dark matter density distribution, a full calcu- 
lation of the gamma-ray emission in the large dimensional parameter space 
of SUSY models is not warranted. Instead, we use the approximate Hill spec- 
trum [38] based on the leading- log approximation (LLA). For MSSM models, 
this approximation is reasonable for neutralino masses below the W mass and 
above the top quark mass. Furthermore, the simplified spectra derived be- 
low vary from the fit obtained in [39] by no more than a factor of a few for 
neutralino masses from 10 GeV up to a few TeV. 

Assuming that the hadronic jets contain only pions and that equal numbers 
of 7r + 's, 7r~'s, and 7r°'s are produced, the neutral pion spectrum is given ap- 
proximately by 



with x n = E n o/m x . Furthermore, the probability per unit energy that a neutral 
pion with energy E n o produces a photon with energy E 1 through the process 
shown in Eq. (5) is 2/E n o. Thus, from Eq. (6) we get for the continuum photon 
spectrum, 




with x 1 = E 1 jm x and y = E n o/m x . The number of continuum photons pro- 
duced per annihilation through this channel with energy greater than a spec- 
ified threshold is found by integrating this spectrum, which gives 



7 



1 dN, 



N cont .(E, > E th ) = J -j^dx, (8) 



x t h 



= l x th 2 - ^T x th + 5v^I + — 1= - ^ (9) 
6 3 by^tt 3 

where x th = E th /m x . 

The 7-ray emission coefficient j is 

j = N,(E > E th )^tLp\r) , (10) 

where iV 7 is the number of photons above the energy threshold as found in Eq. 
(8). The quantity (crv) 1 is the thermally averaged cross section times velocity 
for annihilation into 7-rays, and p(r) is the density profile of the halo. The 
specific intensity of 7-rays, I(R), as a function of projected radius R, is given 
by 

I(R) = j-2 J jdz, (11) 



for which z is the coordinate along the line of sight and z max (R) = y r^ aa , — R 2 
where r max is the radius of the object. The 7-ray flux F from an object at a 
distance d from us is given by 



The radial dependence of j is confined to p 2 {r). We isolate the dependence of 
the flux on the specific halo profile and distance to the LMC by defining K as 

K = 24rf / Ar)*r . (13, 





The 7-ray flux is then given by 

f = K^qh . (14) 

mi 

A. 



For r max we use 3.1 kpc, which is the radial extent of the LMC as observed in 
7-rays by EGRET [40]. We take the distance to the LMC to be 50.1 kpc as 
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found in [23]. We give K in units of (GeV/c 2 ) 2 /cm 5 , which gives the flux F 
in cmT 2 s _1 with (cn>) 7 in cm 3 s _1 and m x in GeV/c 2 . 

For the NFW model K = 3.65 x 10 19 (GeV/c 2 ) 2 /cm 5 , while for the Hayashi 
et al. fits K = (0.84 - 0.85) x 10 19 (GeV/c 2 ) 2 /cm 5 depending on the ex- 
act value of the scale radius. For the isothermal sphere with a core we obtain 
K = 0.43 x 10 19 (GeV/c 2 ) 2 /cm 5 . The more centrally concentrated Moore et al. 
and the Stoehr et al. profiles give slightly higher fluxes, by factors ~ 2 — 4 with 
respect to the NFW case. Clearly the small variations between the fluxes pro- 
duced by these models are overwhelmed by observational uncertainties. Fits 
to all of these profiles give K on the order of 10 19 (GeV/c 2 ) 2 /cm 5 . These pro- 
files behave similarly and consequently give similar fluxes, thus, for simplicity 
we focus on the NFW profile. The mild dependence on the halo profile, and 
more specifically on its central slope, is expected given that we calculate the 
flux from the entire observed extent of the LMC, rather than just from the in- 
nermost region. Thus, the flux depends most strongly on the "normalization" 
of the profile which is set by the mass enclosed in the volume over which we 
integrate. 

EGRET detected a flux of (14.4±4.7) x 10~ 8 photons (E > 100 MeV) cm" 2 s" 1 
from the LMC [41] . Using our fit to the NFW profile and an energy threshold of 
100 MeV, the maximum flux produced by a viable SUSY model is ~ 3.3 x 10~ 9 
photons cm~ 2 s _1 , corresponding to m x ~ 50 GeV and (crv) 1 ~ 2 x 10~ 26 
cm 3 s" 1 . While being consistent with the flux detected by EGRET, even the 
maximum predicted flux is almost two orders of magnitude too low, suggesting 
that the primary source is cosmic rays. Consequently, EGRET's observations 
do not constrain the parameter space. 

EGRET's measurement indicates that cosmic ray induced 7-rays may be an 
additional background component to consider when trying to detect flux from 
neutralino annihilation. Following a method similar to that presented in [42], 
we calculate this cosmic ray induced background. We review very briefly here 
the basic points of the calculation and refer the reader to that study for more 
details. The basic assumption is that cosmic ray acceleration takes place in 
supernovae remnants. Working in the frame of the "leaky box model" for 
cosmic ray propagation, the authors calculate the expected cosmic ray flux in 
the LMC by considering the cosmic ray flux measured for the Milky Way and 
the supernova rates observed for the Milky Way and the LMC. For the 7-ray 
emissivity they include 7-rays from decays of neutral pions produced by proton 
and heavier nuclei interactions, as well as 7-rays produced via bremsstrahlung 
of cosmic ray electrons. 

The calculation of the cosmic ray induced 7-rays in [42] focuses on the flux 
coming from the entire LMC. However, from an observational standpoint it 
is useful to examine the angular dependence of the signal. We calculate the 
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R (kpc) R (kpc) 

Fig. 2. Left panel: Specific intensity as a function of projected distance from the 
center of the LMC for several background contributions and the neutralino signal 
above 1 GeV. The signal is shown for the NFW profile and the isothermal sphere 
with a core (IS). The background components are labeled as follows: G for galactic, 
EG for extragalactic, CRs for cosmic rays. Right panel: Same as left panel but for 
an energy threshold of 50 GeV. The background components in this case are labeled 
as follows: H for hadronic and E for electronic. 

specific intensity as a function of projected distance from the center of the 
LMC. To do this, it is necessary to assume a certain distribution for the HI 
and H2 gas in the LMC, rather than using a mean gas surface density as 
in [42]. According to [26] , the gas in the LMC appears to be distributed 
approximately as a disk with a z-scalelength of about 600 pc and the HI disk 
is about 7.3 kpc in diameter. Assuming an exponential distribution in R as 
well, we equate the 7.3 kpc to 6 radial scalelengths and obtain a scalelength 
of ~ 1.2 kpc, which is close to the scale length of the stellar distribution (see, 
e.g., [28]). Lastly, we normalize the density profile so that when integrated 
over the line of sight coordinate, z, and the projected radius, R, it yields a gas 
mass in agreement with observations. (Note that because the LMC is viewed 
approximately face-on, the rotation axis corresponds approximately to the line 
of sight direction and here we use z to describe both.) 

In Figure 2 we show a comparison of the specific intensity of the different 
backgrounds, as well as that for the NFW profile and the isothermal sphere 
with a core. The left panel assumes an energy threshold of 1 GeV which is 
appropriate for GLAST, while the right panel assumes a 50 GeV threshold 
appropriate for ACTs. For GLAST the relevant backgrounds are the galactic 
and extragalactic diffuse emission. We use the expression found in [43] for 
the isotropic extragalactic background and the expression from [44] evaluated 
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Fig. 3. Angular dependence of the signal-to-noise ratio for a NFW dark matter halo 
profile. For GLAST we use ^4 e // = 10 4 cm 2 and 1 year of observation. For ACTs we 
take A e ff = 10 9 cm 2 and 1 month of observation. 

at the LMC galactic coordinates (I = 280.46°, b = -32.89°) for the galactic 
background. For ACTs the relevant backgrounds are the hadronic and elec- 
tronic cosmic ray shower contributions [44,45] , though it is worth noting that 
ACTs can reject hadronic showers with high efficiency. Clearly, the cosmic ray 
background is not dominant at these energy thresholds. At very low energy 
thresholds (100 MeV - not shown in the Figure) the cosmic ray background 
does become dominant in the very central regions (inner ~ 0.5 kpc). As a 
consistency check, after calculating the specific intensity of the cosmic ray in- 
duced flux above 100 MeV, we integrate it over solid angle and recover the 
~ 10~ 7 photons cm -2 s -1 prediction of [42]. This value is approximately equal 
to the flux measured by EGRET. 

Observationally, the most relevant quantity is the signal-to-noise ratio and its 
angular and energy dependence. In Figure 3 we present the signal-to-noise 
for the flux within an angle 9 , assuming the NFW profile. We show results 
for three different energy thresholds with the neutralino parameters for each 
chosen to optimize the signal. We use the specifications for each instrument as 
indicated. For the noise all the relevant backgrounds, including the cosmic ray 
induced emission, were used. From the figure it is clear that focusing on the 
central regions may help to achieve higher signal-to-noise ratios. As expected, 
for the isothermal sphere with a core (not shown here), the results are less 
optimistic. 

Taking the idealized case of no systematic errors, the detectability condition 
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requires that the signal-to-noise exceed the desired significance s, 

>s, (15) 



Fmin-^-efft 



where F min is the minimum required flux for an s-a detection to be achieved 
by an instrument of effective area A e ff, integration time t, and number of 
background counts for the exposure iV&. The detectability of the neutralino 
induced flux for a particular SUSY model depends on the neutralino mass and 
the annihilation cross section for that model. The neutralino induced flux must 
exceed F min in order to be detectable. This condition can be used to divide 
the m x -(av) y parameter space into detectable and undetectable models for 
each instrument. This is shown in Fig. 4. The points represent SUSY models 
produced using the DarkSUSY package [46]. Models with m x -{av) 1 above a 
given line are accessible to observation by the corresponding instrument, while 
models in the lower region do not yield a detectable flux. 
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Fig. 4. The minimum detectable (crv)^ versus m x for the NFW profile. SUSY models 
above each curve yield a detectable signal for the instrument and observational 
parameters assumed. The dotted line represents observations by upcoming ACTs 
with 50 GeV threshold, 100 hours of observation, and 10 8 cm 2 effective area. The 
dashed line shows detectability by ACTs with an effective area of 10 9 cm 2 which will 
be achievable only at high energy thresholds (~ 1 TeV). The solid line is derived 
for GLAST using the expected flux sensitivity ~ 10~ 10 cm~ 2 s _1 . The dot-dashed 
line is derived for GLAST for 1 year of observation, an effective area of ~ 10 4 cm 2 , 
and an energy threshold of 100 MeV. 

For GLAST, the solid line is derived for the optimistic case of a flux sensi- 
tivity of ~ 1CU 10 photons cm -2 s _1 , which is the expected point source flux 
sensitivity of GLAST for energies above ~ 1 GeV [47]. This flux sensitivity 
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corresponds to a 5-cr detection and 1 year of on-target observation. During 
normal operation, each source will only be in GLAST's field of view about 
20% of the time. Consequently, the actual time required to achieve a detec- 
tion of the same significance is approximately 5 times the necessary on-target 
time. This is within reach since the GLAST science requirement is 5 years of 
operation, while the lifetime goal is 10 years. If GLAST achieves its expected 
sensitivity, then it will be able to detect the neutralino signal for a significant 
portion of the parameter space. This is particularly true in light of the re- 
cently derived limit m x < 500 GeV [4] (vertical arrow in Fig. 4). Additionally, 
a typical value for the total thermally averaged cross section times velocity 
is (av)tot — 2 x 10~ 26 cm 3 s _1 (using Q m h 2 ~ 3 x 10 -27 cm 3 s -1 /(<7v)tot from 
e.g., [2], with Vt m h 2 ~ 0.135 [1]). By making the approximation that all neu- 
tralinos annihilate into pions, and that all neutral pions (roughly one third 
of all pions) annihilate to produce 7s, (en>) 7 is about a factor of 3 smaller 
than (av) tot . Thus, a typical value for (av) 1 is, e.g., ~7x 10~ 27 cm 3 s _1 . If no 
annihilation signal is detected, the corresponding part of the parameter space 
can be ruled out. 

For an energy threshold equal to 100 MeV we take A e ff = 10 4 cm 2 , 1 year 
of on-target observation and Ail LMC ^ 1.13 x 10~ 2 sr for the solid angle 
relevant for the collection of noise (since the LMC will be resolved by GLAST). 
With the sensitivity for this energy threshold only a small portion of the 
parameter space will be detectable (dot-dashed line in Fig. 4). Clearly, the 
best chance for detection is by using moderate energy thresholds, such as 
the 1 GeV case, for which the backgrounds are relatively low but the source 
photons are still numerous. To identify neutralino annihilation as the origin 
of the observed flux, the spectrum and its characteristic features, such as the 
cutoff at E — m x , may be useful. In addition, monochromatic lines produced 
by neutralino annihilation (e.g., the 77 line at E — m x ) can be excellent 
observational signatures for a small region of parameter space where these 
processes are not strongly suppressed (see, e.g., [14]). 

The detectability prospects for existing and upcoming ACTs are less opti- 
mistic. The commonly assumed specifications are A e ff = 10 s cm 2 , E t h = 50 
GeV and 100 hours of observation. For these specifications we find that no 
models will be detectable. A large integration time (~ 1 month) and effective 
area (~ 10 9 cm 2 ) (dotted line in the Figure 4) would improve the chances for 
detection. However, such integration times are fairly long for ACT observa- 
tions and such large effective areas for an energy threshold of ~ 50 GeV are 
beyond the goals of existing and upcoming ACTs. Effective areas of order 10 9 
cm 2 are expected to be achieved by ACTs for energies ~ 1 TeV, but the num- 
ber of 7-rays produced by dark matter annihilation with energies ^ 1 TeV is 
expected to be zero for models where the upper limit on m x < 500 GeV holds. 

Of the upcoming ACTs, the LMC will only be visible to HESS [48] and CAN- 
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GAROO III [49] due to its position in the sky. However, the zenith-angle 
dependence of the effective area of these instruments greatly decreases the 
prospects for detection. The minimum zenith angle of the LMC (at local merid- 
ian crossing) is ~ 46° and 39° for HESS and CANGAROO III, respectively. 
The energy threshold increases and the effective area is drastically reduced for 
these large zenith angles, requiring even longer integration times. 

In this section we focused our flux calculations on the NF W fit and also dis- 
cussed the less favorable case of the isothermal sphere with a core. In addition, 
we estimated the change in flux for other profiles. Of the range of profiles we 
considered, the fluxes produced vary by less than an order of magnitude, and 
hence the detectable cross sections and neutralino masses do not differ greatly 
from those of the NFW profile. This mild dependence on the detailed profile 
arises from our integrating to a radius R which includes most of the LMC as 
it appears in 7-ray observations. We find that the cosmic ray induced back- 
ground is largely subdominant compared with other background sources, and 
that for the signal to exceed all backgrounds, observations of the central region 
of the LMC may be the most promising. 



4 The synchrotron emission 

Neutralino annihilation not only generates neutral pions, but also a compara- 
ble number of charged pions. The charged pions decay as 

7r + — > /U + z/ M and 7r~ — > /i~z? M . (16) 

Muons then decay via 

fi + -> e + u fl u e and fi~ -> e"z/^z/ e . (17) 

In the presence of magnetic fields, the electrons and positrons produced in 
this way generate synchrotron radiation. The synchrotron emission from neu- 
tralino annihilation has been considered for a variety of objects including the 
galactic center [6], dark matter clumps [17], Draco [10], and halos of galaxy 
clusters [22]. Here we follow previous studies and calculate the synchrotron 
emission due to neutralino annihilation in the LMC. 

The radio flux from the LMC has been measured in frequencies from 19.7 
to 8550 MHz. The electron energies with maximum synchrotron emission in 
this range are relatively low (~ 20 GeV for the most energetic 8550 MHz 
synchrotron photons for a 5 /xG magnetic field). Thus, we can consider only the 
dominant source of electrons and positrons for these energies: the 7r ± decays. 
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This is valid for gaugino-like neutralinos (see [22] and references therein). 
We use the analog of Eq. (6) (substituting tt° with charged pions which we 
denote by ir) for the charged pion spectrum, dN n /dE w , produced by quark 
fragmentation. 

The number of electrons and positrons produced per annihilation per energy 
interval dN e /dE e is then given by 

dN m f a fdN dN^dN^ 

1 1 ' 1 " dE w dE^ (18) 



dE e J J dE n dE^ dE e 



where a = (m^/m^) , 

^ = 1 < (19) 
<i£„ E, ml -ml' { ' 



and 



d/Vg 
dE, 



E„ 



3 Ee 

2 \E„ 



+ 



2 (Ee 

3 



(20) 



Eqs. (19) and (20) give the decay products from charged pion and muon 
decays, respectively. For simplicity, we adopt the notation n for n^ 1 and e for 

e ± / 

The source spectrum q e which is the number of electrons per unit time, volume, 
and energy can be written as 



with (<jv ) in this section denoting the thermally averaged cross section times 
velocity for annihilation into charged pions. For electron energies relevant to 
the frequency range that we use, the source spectrum can be altered by syn- 
chrotron losses and inverse Compton scattering (ICS) (see, e.g., [50]). The 
final spectrum dn e /dE e including losses is 



dn P 

q e r (22) 



dE e 
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where r is the average lifetime of an electron with energy E e , 

T ~ dEjdt ■ (23) 



The dominant loss process and, thus, the shortest lifetime, is used in Eq. (22). 

The total synchrotron power emitted from an electron of energy E e , in the 
presence of a magnetic field = \B\ in microgauss, is [51] 

= 1.6 x KT 15 ^ 2 erg/s . (24) 

(XL 



Combining Eqs. (23) and (24) we find r syn , the time scale for energy losses 
due to synchrotron emission. 

In order to estimate r syn one needs the local magnetic field B^. There are nu- 
merous estimates in the literature of the component of the magnetic field along 
the line of sight B^^ os . Rotation measure analyses of pulsars in the LMC indi- 
cate a B^ i os ranging from ~ 0.4 to 5 fiG [52,53]. Similar values {B^^ os ~ 2 fiG) 
were obtained from the polarized radio continuum emission [54] . Estimates of 
the total magnetic field, rather than the component of the magnetic field along 
the line of sight, are more indirect and typically involve more assumptions. 
Several studies seem to converge to ~ 6 fiG [55,56], while the largest value es- 
timated for the total magnetic field is ~ 18.4 fiG [57]. In what follows we adopt 
the value of ~ 5 fiG for the total magnetic field, unless otherwise specified. 

For magnetic fields of 5 fiG, synchrotron losses dominate over ICS processes 
off the cosmic microwave background, the optical photons, and the produced 
synchrotron photons themselves. Therefore, the relevant timescale in Eq. (22) 

IS Tgy n . 

The synchrotron emission coefficient j u (number of photons of frequency v 
emitted per unit time, volume, and frequency) is 

. _ dn e dE e dE syQ erg 

3u ~ dE e dv dt cm 3 sHz ' [ ' 



To calculate the factor dE e /du appearing in the expression for the emission 
coefficient, we use the relation between E e and the frequency where maximum 
synchrotron emission occurs for this electron energy. The frequency of max- 
imum synchrotron emission is equal to ~ 0.29z/ c , with the cutoff frequency 
v c = 3/27 3 z/ B where 7 is the Lorentz factor and v B = (eB sin a/2-K^mc) is the 
relativistic gyrofrequency. Taking the power weighted average over the pitch 
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angle a, the frequency of maximum emission is 

E e \ 2 



^ = 4.1^(^) MHz. (26) 



Using Eqs. (18), (21), (22), (23), (24), and (26) to calculate j u from Eq. (25), 
and substituting the latter in Eq. (12) we calculate the synchrotron flux. For 
the integration we use the same R that we use in the 7-rays since the extent 
of the LMC in the two frequency ranges is similar [40]. 

The synchrotron flux F syn in Jy (1 Jy=lCT 23 ergs cm~ 2 s _1 HzT 1 ) is 



where K is given by Eq. (13) and is in units of (GeV/c 2 ) 2 cm~ 5 . Again, the 
neutralino mass m x is given in GeV and (av) in cm 3 s _1 . The frequency 
dependence of the flux is given by 



— - oc 2.48 - 1.22aT - 5 + 0.22a; 2 + O.lOa;- 1 - 5 - 1.54a; - 5 - 0.04a; 3 , (28) 
dE e 

where x = E e /m x c 2 ~ (4.94 x 10~ 4 -B~ a5 z/ a5 /m x ), with v in Hz and m x in 
GeV. The dependence on v varies with the electron energy and the neutralino 
mass and goes as v~ - 75 at sufficiently small energies. 

A compilation of radio observations of the LMC from 1959 to 1991 in the 
frequency range (19.7 — 8550) MHz is presented along with the corresponding 
references in [56]. As can be seen in Fig. 5, the observed LMC spectrum 
in the (19.7 — 8550) MHz range has a steep contribution from synchrotron 
radiation which dominates at lower frequencies until around 2 x 10 3 MHz, at 
which point the thermal emission becomes dominant. A least squares fit of 
the form S u = S 0tt hV~ 0A + S 0jnt hV~ anth to the data yields a total flux density 
at 1 GHz of 610 ± 10 Jy, a non-thermal spectrum with a nt h = 0.70 ± 0.06, 
and a fraction of thermal emission at 1 GHz of (45 ± 10)% [56]. Using this 
information we find S ,th — 2.2 x 10 3 and S , n th — 6.7 x 10 8 . The data points 
and the decomposition into thermal and non-thermal emission are shown in 
Fig. 5. We also plot F syn for m x = 50 GeV and (av) = 2.2 x 10~ 26 cm 3 s~\ for 
the NFW dark matter profile. The thick solid line corresponds to = 5/iG; 
the lighter line represents the results for the maximum total magnetic field 
estimated in literature (18.4 fxG [57]). 

In the more realistic case of B^ = 5 fj,G, neutralino-induced synchrotron emis- 
sion clearly may be part of the observed flux, but is lower than the total 
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Fig. 5. The flux in Jy as a function of frequency in the range from 19.7 to 8550 
MHz. The open triangles correspond to LMC data taken over the period from 1959 
to 1991; the dashed line is a best fit to the data assuming a decomposition into 
thermal and non-thermal emission. The dotted line corresponds to the non-thermal 
contribution, and the dot-dashed line to the thermal contribution. The flux from 
neutralino annihilation for m x = 50 GeV, = 5 //G, and (av) ~ 2x 10 -26 cm 3 s _1 , 
is given by the thick solid line. The remaining line corresponds to = 18.4 fiG, 
which is the maximum published estimate for the total magnetic field (see text for 
details). The B^ ~ 5 p»G case is the most realistic value, even though the results in 
these frequencies do not change dramatically even when a 18.4 [iG field is used. 

observed flux, especially at higher frequencies. In the low frequency region, 
the signal increases enough to exceed the thermal contribution to the obser- 
vation. Given that the neutralino signal (oc z/~ a75 ) is almost parallel to the 
non-thermal component (oc z/ -aro ), even at reasonably low radio frequencies 
the neutralino flux is sub-dominant. A mildly constraining limit on neutralino 
cross-sections can be reached by requiring that F syn be less than or equal to 
the non-thermal contribution in S u at low frequencies. The frequency depen- 
dence essentially drops out and we find (erf) < 1.24 x 10~ 9 m^ 5 B^ 25 j K . Thus, 
assuming the NFW fit and m x = 50 GeV, any (av) lower than ~ 2 x 10 -26 cm 3 
s _1 gives a flux at low frequencies consistent with observation. The allowed 
(av) become even larger at larger neutralino masses. Thus, the observations 
do not put severe constraints on the allowed cross sections. Although we have 
focused on the NFW profile, the synchrotron fluxes vary slightly depending 
on the halo profile used, as discussed in §111. 

The best hope for distinguishing the synchrotron emission generated by cosmic 
rays from that produced by neutralino annihilation is to use the fact that the 
density profile of cosmic rays differs from the dark matter density profile. 
While the dark matter halo extends significantly beyond the disk, the density 
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of cosmic rays is expected to rapidly decrease at large radii. If low radio 
frequency observations of the LMC are made with high angular resolution, the 
change from cosmic ray to neutralino dominance should be apparent as one 
moves away from the LMC disk. Currently there is no telescope that could 
carry out observations of the LMC at frequencies less than about 20 MHz. 
Such low frequencies are extremely difficult to observe from the ground due to 
ionospheric absorption and scattering. One promising ground-based project 
is the Low Frequency Array (LOFAR) [58] which should reach frequencies 
down to ~ 10 MHz and a flux sensitivity of a few mJy in 1 hour. However, 
the site of the instrument has not yet been decided, so the LMC may not 
necessarily be observable. To reach even lower frequencies, where ionospheric 
absorption is very intense, space-based instruments are required. The proposed 
Astronomical Low Frequency Array (ALFA) [59] should reach down to ~ 0.3 
MHz. 



5 Discussion and conclusions 

We have studied the 7-ray and synchrotron emission from neutralino anni- 
hilation in the LMC. We modeled the dark matter profile of the LMC by 
fitting modern profiles found in N-body simulations to rotation velocity data. 
Although we focused on the NFW best fit density profile when deriving our 7- 
ray and synchrotron fluxes, we also considered a less concentrated dark matter 
profile (the isothermal sphere with a core) and a more concentrated one (the 
Moore et al. profile). This range of profiles changes our results by less than an 
order of magnitude for the same rotation curve. If one takes into account the 
expected compression of the dark matter due to baryons, which would cause 
even a shallow, core-dominated dark matter distribution to become cuspy, it 
is reasonable to assume that the actual fluxes may be several orders of mag- 
nitude greater. Prada et al. [33] found that baryonic compression can enhance 
the central density by more than an order of magnitude. 

The LMC rotation curve data has been recently re-analyzed by van der Marel 
et al. [23] who claim that the plateau value for the rotation velocity is not 
~ (60-70) km/s as found in most studies, but ~ 50 km/s. Their conclusion 
is based on a re-analysis of carbon star data including corrections for the 
transverse motion of the LMC, as well as for nutation and precession effects. 
Their resulting rotation curve has large errors and is sparse at small radii, 
but a fit to their rotation curve results in K up to ~ 1 order of magnitude 
less than those we found using the initial data. The quality of the fits is much 
poorer given the sparseness at small radii. 

An additional uncertainty is introduced by the fact that only the line-of-sight 
component of the velocity can be observed. The assumptions necessary to re- 
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cover the total velocity from the line-of-sight component may result in an over- 
or underestimation of the actual velocity, thereby improving or worsening the 
prospects for detection. 

The predicted cosmic ray induced 7-ray and synchrotron emission appears to 
be of similar magnitude to the observed emission. A high angular resolution 
study of both neutralino and cosmic ray induced emission from the LMC may 
be key to utilizing the closest known dark matter "clump" as a laboratory for 
testing neutralino parameters. Smaller solid angles help reduce the collection of 
noise, and, most importantly, a study of the spatial dependence of the emission 
can help GLAST and radio observations distinguish between neutralino and 
cosmic ray produced signals. 

Summarizing, we find that neither the EGRET measurement of 7-rays above 
100 MeV from the LMC nor the low frequency synchrotron observations signif- 
icantly constrain the SUSY parameter space. The maximum 7-ray flux for the 
considered profiles is almost ~ 2 orders of magnitude less than the observed 
flux, whereas all models with (cry) ^ 2 x 10~ 26 cm 3 s _1 are consistent with 
the observed synchrotron flux. Existing and upcoming ACTs will not be able 
to probe the SUSY parameter space, unless the neutralino profiles are more 
centrally concentrated or observers use longer integration times and larger 
effective areas larger than currently planned. Since the sensitivity becomes 
better for a larger exposure (effective area x integration time), the problem 
of a small effective area may be compensated for by larger integration times. 
However, these integration times may not be realistic. Added to these diffi- 
culties, the location of the LMC makes it difficult to observe even for some 
instruments in the southern hemisphere. 

GLAST on the other hand, being a satellite, will not find the LMC location 
problematic. If systematic errors are small, GLAST will be able to probe a 
significant part of the allowed parameter space, especially if the recent limit 
m x < 500 GeV is taken into account. The comparison of the 7-ray annihilation 
signal to the EGRET measurement implies that there is another dominant 
source of this flux, most likely cosmic rays. The synchrotron study finds this 
result as well. Thus, the spectral features of the neutralino annihilation signal, 
in particular the shape and cutoff of the spectrum at m x , will be essential to 
disentangle the neutralino signal. 

The synchrotron signal is at best a factor of ~ 3 less than the observed signal, 
with this discrepancy increasing with increasing frequency. The flux values for 
the frequency range we use vary from ~ 1 to ~ 1000 Jy, easily detectable 
levels. The difficulty will again be in disentangling the neutralino induced 
signal from the total flux. At low frequencies, the frequency dependence of the 
calculated flux is almost identical to that of the observed spectrum. The best 
hope for observing the synchrotron emission due to neutralino annihilation is 
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a high angular resolution study at low frequencies which may be possible with 
proposed observatories such as LOFAR or ALFA. 
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